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Abstract 



We have studied the superfluid - Mott insulator transition in a BEC con- 
fined in an optical lattice and driven by an additional external electromagnetic 
field, which excites another internal state. It is shown that due to resonant 
strong driving the critical values of the transition parameter depend on ex- 
ternal field parameters (the Rabi frequency and the detuning from resonance) 
and can be tuned to a given value. It is also shown that for driven system 
there are two optical potential depths at a phase transition, which are always 
more than without driving. 

* Electronic address: ggenkin@argento.bu.edu. 



1 



Typeset using REVT^X 



The experimental observation of Bose-Einstein condensation in a dilute gas of ultracold 
trapped atoms [ 1-3 ] has generated much interest in manipulating such coherent matter 
by external electromagnetic fields. Recently, atoms have been confined in optical potentials 
created by standing light waves [ 4, 5 ], with the wavelength of the optical potential much 
smaller than the dimensions of the trap. Considerable work has been done, in the last decade, 
to determine the ground state phase diagram of correlated bosons on such a lattice. The 
bosons can undergo a transition from superfluid to insulator. Strong competition between 
potential and kinetic energies begets a quantum phase transition. The on - site repulsion 
can produce a Mott insulating phase with a quantum phase transition to a superfluid as 
the interaction strength weakens. A seminal experiment by Greiner and collaborators [5] 
demonstrated a quantum phase transition in a Bose - Einstein condensate ( BEC ) from 
a superfluid ( SF ) state into a Mott insulator ( MI ) state, by varying the lattice laser 
intensity, as proposed theoretically in [ 6 ]. The quantum phase transitions have attracted 
much interest in recent years. These transitions are accomplished by changing not the 
temperature, but some parameter in the Hamiltonian of the system; for a BEC in an optical 
lattice this parameter is the lattice laser intensity, for solid state systems, for example, the 
magnetic field in a quantum - Hall system, or the charging energy in Josephson - junction 
arrays. 

In this paper, we propose a manipulation of the quantum phase transition of a BEC 
in an optical lattice by an additional external electromagnetic field. This basic idea of my 
paper is using the external field to drive transitions between internal states for atoms in 
an optical lattice. This driving alter the balance between self - interactons of the atoms in 
lattice site and hopping of the atoms between the sites, and thereby modify the superfluid 
- Mott insulator phase transition. 

An applied resonant field excites another internal state of the BEC, which is why the 
populations of two coupled states are varied. Without driving a quantum phase transition 
in an one - component BEC is a single - channel problem, but with resonant driving it 
becomes a two - channel problem, in which every channel has competition between kinetic 
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and potential energies. The populations of the coupled states are varied. As a result, the 
effective on - site interaction and the effective tunneling in these two channels, which depend 
on the population, are also changed. However, due to this varying population the variations 
of the on - site interaction and the tunneling are different; therefore, the balance between 
kinetic and potential energies, which determines a phase transition, is altered. As a result, 
we have a modification of the quantum phase transition by the applied resonant field, which 
is determined by the two channels. It is shown that for each channel, the critical value of the 
ratio of the interaction strength and the tunneling depends on the population, which for the 
strong driving is the Rabi oscillations. Therefore, in the regime of the time - dependent Rabi 
oscillations the MI - SF and the reverse SF - MI transitions occur during a period of the 
Rabi frequency. The frequency of this oscillating superfluidity is determined by the driving 
strength. In the steady state driving regime the critical value of the ratio of the parameters 
is always increased, and two optical potential depths at the phase transition for the driven 
system are always higher than without driving. Therefore, this effect opens the possibility 
creating a superfluid state from a initially Mott insulator state. The critical values of the 
driven system depend on the parameters of the external field ( the Rabi frequency and the 
detuning from resonance ); in particular, the quantum phase transition can be induced by 
varying only the detuning with the fixed field amplitude. We can tune the critical values of 
the quantum phase transition to a given value. Our proposal of manipulating a quantum 
phase transition can be applied to solid state systems. For example, in a quantum - Hall 
system the resonant driving may be a novel way of modification a quantum phase transition 
besides using the magnetic field which controls the transition between different quantized 
Hall phases. 

We consider a two - component Bose - Einstein condensate in an optical lattice in the 
presence of an additional external electromagnetic field which couples the two components 
(the different internal states). The energy eigenstates are Bloch states, and an appropriate 
superposition of Bloch states yields a set of Wannier functions which are localized on the 
individual lattice sites. An external field couples the two states a and b, and an atom is always 
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in the superposition of these states. The system is driven by a two - photon Raman pulse, 
the strength of which we characterize by the Rabi frequency Qr. We will consider as the two 
distinct internal states a and b coupled by the external time - dependent field two different 
vibrational states in the Wannier basis, where the initial state a is the lowest vibrational 
state of the harmonic oscillator with n = and the state b is the excited vibrational state 
with n^O. Depending on the relative strength of the atom - atom interaction, the system 
may lie in the Rabi regime or in the Josephson regime. The equations of motion of a 
strongly ( the Rabi regime ) driven [ 7, 8 ] two - component BEC resemble the Bloch 
equations describing a driven two - level system if the external driving field MIr is much 
larger than the difference in the two - body interaction ( the " charging" [ 9 ] energy E c 
which is proportional to U a + Uf, — 2U a b ).For the two - body interaction, we use a contact 
potential, V;(x — x') = [/^(x — x'), here {1 = a, b}; the coupling constants, Ui, are expressed 
in terms of the s -wave scattering lengths by U\ = ^ n m ai , where a a and a b are the scattering 
lengths for collisions between two atoms in states a and b, respectively, while a a b is the 
scattering length for collisions between atoms in different internal states. The driving field 
is a spatially uniform field of the two - photon pulse of Raman type with the wavelengths 
much larger than the dimensions of a BEC. The external resonant field excites the higher 
internal state, while, the energies involved in the one channel dynamics to be small compared 
to the excitation energies to the higher internal state. In this case of the strong driving we 
can solve the problem in two steps. The first is taking account only of the strong external 
field Qr, in the result of this step we have the Bloch equations solution tpi where the functions 
ipi are determined by the parameters of the external field ( the Rabi frequency Qr and the 
detuning 5 from resonance ). This Rabi solution of the Bloch equations obey the condition 
|<Aj| 2 + I<a| 2 = 1- And, the wave function of a two - level system is 

= ¥>a*o + ¥>6*6- (1) 

In this two - state Hilbert space the system is described by the effective Hamiltonian H e ff. 
And, a second step is taking account of the remaining terms of the Hamiltonian. As a result, 
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using Eq.(l) we can write the effective Hamiltonian in the second - quantization notations 
for the boson field operators 

H eff = J ^[*t( X) t)(-^ + y (xH^ 

J d 3 xd 3 x'[S>l(x, *)*1(X', t)K(x - x')* a (x', *)* (X, t)<pi+ 
tfJ( X , t)tf J(x', t) V 6 (x-x')* 6 (x', *)* 6 (X, t)4+2¥ a ( X , t)*t(x!, t)K 6 (x-x')*6(x', i)* a ( X , t)¥>M, 

(2) 

where Vo(x) is the optical lattice potential, usually, the optical lattice potential has the 
form Vo(x) = Y%=i Vj sin 2 kxj with wave vectors k — ^ and A the wavelength of the laser 
beams. The atomic field operators have been written in the rotating wave approximation 
and, consequently, in Eq.(2) there appears a detuning, 5 = u — u d , where Tiujq is the 
energy difference between the two internal states, and Ud = ujl — us is the frequency of 
the Raman drive. Following Ref. [6] the boson field operators can be represent by an 
expansion in the Wannier basis, which is a sum over the lattice i, ^/ a ( x ) — Si a i w o( x — 
Xj),^6( x ) = J2ibiW n (yL — Xj), where the operators ai and bi are the operators on lattice 
site i; the annihilation and creation operators a* and af ( bi and bf ) obey the canonical 
commutation relations [a*, Oj] = Sij([bi,bj] = 

As a result, in the tight - binding approximation a two - band Bose - Hubbard Hamilto- 
nian of a strongly resonant driven BEC is 

H B h = -*«» Y a \ a W 2 a - hb Y h \ h wl 

<i,j> <i,j> 

+^7p E a Ui(ahi - l)Va + ¥jr Y bibi(b\bi - l)<p\ + C/ab I] a\aib\bi(pl(t)(pl+ 

i i i 

a\a^ 2 a + Y W) + ^ E (3) 

i i i 

where the sum in the two terms runs over all nearest neighbors denoted by < i,j >. The term 
involving the chemical potential \i is added as the calculation is in the grand - canonical 
ensemble at temperature T = 0. The first and second sum measure the boson kinetic 
energy [10], and the parameter tu is the hopping matrix element for channel / between 
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neighboring lattice sites the parameters U aa , U bb , U ab characterize the strength of the on 
- site interaction of two atoms on the lattice site i. These parameters are defined as 

r h 2 

taa = J <FxWq(x - X i)[-^ V2 + V (x)]w (x - Xj), 

r h 2 

t bb = J d 3 xw* n (x - x i)[~7^ v2 + v o(x)]w n (x - Xj), 

Uaa = U S a J d 3 x\w (*)\\U bb = U S b J d 3 x\ Wn (x)\\U ab = J d 3 X \w (x) |>„(x) | 2 . (4) 

The parameters too and [Too of the one - band Bose - Hubbard model depend on the potential 
depth V , and they were calculated analytically in Ref. [11] in the tight - binding limit in terms 
of the microscopic parameters of the atoms in the optical lattice using the approximation 
[11] the Wannier states by the harmonic oscillator wave function u> n (x — Xj) of number n. 

In order to describe the zero - temperature quantum phase transition from the superfluid 
to the Mott insulator phase analytically, it is neccesary to make some appropriate mean - field 
approximation in using the Hamiltonian ( Eq.(3)). A mean - field theory, developed in [11], 
constructs a consistent mean - field theory by substituting a\a,j = {a\)dj + a\{aj) — (a\)(aj), 
where (a\) = (oj) is the expectation value on site i. Following this theory [11], which 
exploited the usual Landau procedure for a second - order phase transition, we obtain the 
critical value of the ratio of the on - site interaction and hopping parameters 

(^) (cr V 2 = ma + i + v / (2^ + i) 2 -i), (5) 

where the number of nearest neighbors 2d, d is the number of dimensions, g is the integer 
number of particles at each site i. In these two equations the functions ipi coupled by the 
normalization condition. For a single channel ( one - component model ) ipo = 1 and these 
equations are reduced to the well - known result for the critical value [11]. Note that the 
chemical potential /z c at the critical point is calculated also by the cross - interaction U ab 
between the channels. 

Let us first consider the time - dependent regime of the Rabi oscilations for zero detuning, 
5 = 0, in this case we have the time - dependent functions (pi(t). In this regime the Eqs.(5) 
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are valid if the Rabi oscillations period is larger than the nonequilibrium dynamical times 
(the characteristic time scales are the time of establishing of the SF state, the restoration 
time of the phase coherence or the dephasing time). The nonequilibrium dynamic time for 
trapped bosonic atoms in an optical lattice potential was considered in Ref. [12] and was 
shown that the characteristic time of the dynamical restoration of the phase coherence in 
the domain of a quantum phase transition is the integer of the order of 10 multiplied by a 
Josephson time h/t. The initial phase of the BEC is MI phase for V >- Vq t , where V^' is 
the optical potential depth at the SF - MI transition without driving. The functions ipi(t) 
determine the population pi of the states /, which oscillate as cos 2 Q R t for channel a and 
sin 2 £l R t for channel b. Accordingly, during the period 2ir/il R of the Rabi oscillations the 
transition from MI to SF phases and the reverse transition occur. Really, for the channel a 
at the point 61 there is the MI - SF transition, at (tt — ©1) - the opposite SF - MI transition. 
For times of the Rabi oscillations period from t — until t\ = Q\/Vt R and from (tt — ©i) 
until 7r the channel a is in the MI phase, from 61 until (n — ©i) there is the SF phase. 
Correspondingly, for the channel b with the opposite time - dependence of population at @ 2 
there are the SF - MI transition, at (it — 2 ) the MI - SF transition. In the channel b, until 
t 2 = ©2/^R and from (n — 2 ) until ir is in the SF phase, from 6 2 until (tt — 6 2 ), there is 
the MI phase. Thus, both channels have the superfluid - Mott insulator transition twice in 
each half - Rabi period. The frequency of this oscillating superfluidity is determined by the 
driving strength. The time - averaged population N dr of SF phase of a driven BEC, which 
without driving is in the MI phase, is 

N dr (V ) tt + 2[e 2 (y ) - e 1 (V )} - [sm2© 1 (y ) + sin2© 2 (y )] 



iV 2tt 
where N is the concentration of a nondriven BEC in the SF phase for V -< V cr . Here 



(6) 



@i(^o) — arccos , 



AC/ o(Vb) n 

, B 2 (Vo) = arcsm 



nn 



AC/oo (V ) AC/00 



(7) 



where AU 00 (V ) = ^n,At 00 (V ) = the ratios, which are calculated in Eqs.(10) 

below for n = 1, are AC/ 7 = ^, At® = ^. If lg °(Jg %E > 1 then & 2 (V ) = vr/2. 
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In the steady state driving regime every atom is in the stationary superposition of two 
states, which populations are determined by the Rabi frequency and the detuning from 
resonance. We have two critical values of the transition parameters of the driven system 
{U n /t n )t\ { 1 = a, b }. The dependence Uu(V ),t u (V ) for every channel determines two 
optical potential depths and V^jT at a phase transition, correspondingly. The critical 

value of parameter of a phase transition for the driven system (Uu/tu)^ is always higher 
than without driving (Uu/tu)^. For the channel a 

§ 2 ^C = 1 -5[ 1 + (cr) a i- 1 . (8) 

Therefore, there are an interval of the MI phase, in which the driving switches the MI phase 
to the SF phase. In this band from (t/ooAoo)^ to 2(£y 00 /t 00 )Q Cr * ) , by varying the parameter 
5/Qr of the external driving, we can tune the critical value to a given value. With increasing 
the depth of the optical potential, the on - site interaction increases and the tunneling matrix 
element is reduced. Accodingly, with increasing the critical value of the transition parameter, 
the optical potential depth V^r) a t a phase transition always increases. 
For the channel b the critical value is always more than for the channel a 

(UooMdr' 

We obtain for n — 1 and in the tight - binding limit, A ^> ft, the ratio 

£H*<=4>- (10) 



Here (3 = ^^mv~- * s ^ ne s * ze °^ the g roun d - state oscillator wave function determined by 
the atomic mass m and the oscillation frequency in the wells ( see, for example [ 6 ] ) 



Vj = ^jAEjiVj /h, Er is the recoil energy, H 6 (x) is Hermite polynomial. In general, the ratio 
Unn/Uoo contains sum over the Hermite polynomials H n (0) till 4n+2, and t nn /t 00 ~ (\/2(3) 2n , 
For all n there are the inequalities t nn ^> t 00 and U nn ^> Uqo, this increase t nn and U nn is due 
to oscillations of the wave function for n ^ 0. For example, tn/too ~ 64 and Uu/Uqq ~ 45 
for sodium [13] where E R /h ~ 2n x lO^Hz for a red detuned laser with A ~ 10 3 nm and 



V ~ 16E R . Although, the on - site interaction on the higher channel U nn increases with 
the number n, the tunneling matrix element t nn increases faster, and we can say that the 
higher channel n is " more superfluid". The ratio U nn /t nn increases also with increasing the 
optical potential depth Vb, and because there are the inequalities U nn /t nn -< t/ooAoo for the 
given Vo, then Eq.(9) yields a relationship V^™. >- 

In summary, the resonant strong driving by a spatially uniform field modify the SF - MI 
transition in a BEC in an optical lattice by the mechanism of the turning on of the second 
channel. Therefore, there are the modifications of the main channel and the second channel. 
As a result, the transition on each channel is controlled by the external field. In the time 
- dependent driving regime the initial MI phase switches over to the SF phase and the SF 
phase to the MI phase during the Rabi oscillations period. The frequency of the oscillating 
superfluidity is determined by the driving strength. In the steady state driving regime, the 
critical values of the ratio of the on - site interaction and the atomic tunneling rate depend on 
the external field parameters ( the strength and the detuning from resonance). Two critical 
values of the driven system and, accordingly, two optical potential depths at the phase 
transition are always more than without driving. The proposed mechanism of manipulating 
and controlling Bose - Einstein condensates makes it possible to tune the critical value of 
parameters of the quantum phase transition to a given value. 

I thank H. E. Stanley for encouragement. 
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